Effect of Semicore Orbitals on the Electronic Band Gaps of Si, Ge, and GaAs within 

the GW Approximation 
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We study the effect of semicore states on the self-energy corrections and electronic energy gaps 
of silicon, germanium and GaAs. Self-energy effects are computed within the GW approach, and 
electronic states are expanded in a plane-wave basis. For these materials, we generate ab initio pseu- 
dopotentials treating as valence states the outermost two shells of atomic orbitals, rather than only 
the outermost valence shell as in traditional pseudopotential calculations. The resulting direct and 
indirect energy gaps are compared with experimental measurements and with previous calculations 
based on pseudopotential and "all-electron" approaches. Our results show that, contrary to recent 
claims, self-energy effects due to semicore states on the band gaps can be well accounted for in the 
standard valence-only pseudopotential formalism. 

PACS numbers; 71.15.-m, 71.15.Ap, 71.15.Dx, 71.10.-w 
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I. INTRODUCTION 



Since the early applications of the GW method to 
real materials (see Refs. 1-3 and References therein), 
the pseudopotential plane-wave approach has been the 
method of choice due to its accuracy and technical sim- 
plicity. Recent advances in LAPW and LMTO method- 
ologies have allowed the implementation of "all-electron" 
applications of the GW method4i^iSiL& One common fea- 
ture of such calculations, using standard level of approx- 
imation for the self-energy, is an underestimation of the 
electronic energy gap compared to experimental measure- 
ments, whereas pseudopotential-based calculations show 
very good agreement with experiment iSiS To explain 
this inconsistency, it was proposed that the pseudopo- 
tential approach does not correctly describe the effect of 
core orbitals in the self-energy corrections to the energy 
gaps, resulting in overestimated corrections^iS 

It is thus desirable to elucidate the effect of core or- 
bitals in the quasiparticle band structure, and the pre- 
ferred procedure is to perform a well converged all- 
electron calculation and compare its results with simi- 
larly converged pseudopotential-based calculations. Ob- 
viously, numerical precision should not be neglected. In 
this work, we explicitly include semicore orbitals in the 
pseudopotential plane-wave approach and calculate the 
quasiparticle energy gap for three semiconductors of tech- 
nical importance: Si, Ge and GaAs. The underlying de- 
scription of the ground-state electronic structure is based 
on Density Functional Theory in the Local Density Ap- 
proximation (DFT/LDA).^*y^ Throughout this work, we 
are careful to converge all results systematically, and the 
final results are compared to previous pseudopotential 
results and to recent all-electron calculations. The pa- 
per is organized as follows: we outline the theoretical 
method in section II. Results are presented in section III 
and discussed in section IV. 
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FIG. 1: Convergence of the minimum gap in Si as function of 
the number of unoccupied bands included in the calculation of 
the Green's function, Eq. Q (open triangles). The solid line 
is a guide to the eye. Results obtained by Ku and Eguilua- 
are shown as black diamonds. 



II. THEORETICAL METHOD 

The pseudopotential formalism has two advantages 
that make it convenient for practical ab initio calcula- 
tions. First, degrees of freedom due to core electrons are 
removed from the system, resulting in a description that 
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contains only valence electrons. Second, valence elec- 
tronic pseudowavefunctions are smooth in the vicinity of 
atomic sites as well as in the interstitial region. Such 
valence wavefunctions can be expanded to convergence 
easily in a relatively small basis set such as a plane-wave 
basis. On the other hand, one possible problem with this 
formalism is that using pseudopotentials, instead of the 
true electron-ion potential, may not fully describe all ef- 
fects produced by interactions between the valence and 
core electrons. In order to address this issue in GW calcu- 
lations, we include in the present study electrons from the 
outermost valence shell as well as those from the second 
outermost ("semicore") atomic shell as active "valence" 
electrons in the pseudopotential formalism. Only inter- 
actions between these "valence" electrons and electrons 
from the deeper core shells are described by pseudopoten- 
tials. Taking Si as an example, its "core" now contains 
only Is electrons. Electrons from the 2s,2p,3s,3p, and Sd 
shells are all treated on equal footing in the subsequent 
calculations. 

Apart from the atomic configuration, we follow 
the standard prescription for generating ab initio 
pseudopotentialsiiSii^ Since the deeper core electrons 
have extremely large binding energy {e.g., 130 Ry for 
the Is electrons in atomic Si), their interaction with va- 
lence electrons is expected to be much weaker than the 
already small interaction between the outermost valence 
and the semicore electronsii^ 

The Kohn-Sham DFT formalism within the LDA is 
used to solve for the ground state electronic structure 
and to provide a starting point for the calculation of the 
electron self-energy. We follow closely the GW method as 
developed by Hybertsen and Louiei In this method, the 
self-energy is given by the standard GW approximation, 

S(r,r';£;) = z J ^e'^'^''^^ Go{r,r'; E-E')Wo{r,r'; E'), 

(1) 

where Go is the one-electron Green's function, calculated 
from LDA energy eigenvalues and eigenstates {S = 0+ for 
occupied states, S = 0~ for unoccupied states). 

The screened Coulomb interaction Wq is calculated 
within the Random Phase Approximation (RPA) as 
Wo = [1 - VPo]~^ V, with V being the bare Coulomb 
interaction and the polarizability Pq = —iGoGo.^^ This 
is the commonly employed level of approximation for GW 
calculations {i.e. neglecting self-consistency and vertex 
corrections)"'?, and we specifically compare results from 
both pseudopotential and all-electron calculations at this 
particular level of approximation to help untangle the ef- 
fects of core states. 

The convolution integral in Eq. is performed us- 
ing the generalized plasmon-pole (GPP) model, ^ which 



TABLE L Atomic parameters used to generate semicore pseu- 
dopotentials. Pseudowavefunctions were defined from the or- 
bital with lowest principal quantum number at each angular 
momentum channel. Cutoff radii are given in units of Bohr 
radius. 



Reference Configuration rcut local channel 

s p d 

Si 2s^2p^3s'^3p^'^-'3d°°-' 0.40 0.35 0.40 d 

Ge 3s^3p^3d^°4s'^-ip^ -'4d^-^ 0.50 0.50 0.50 s 

Ga 3s^3p^3d^°4s24pMd'^ 0.50 0.50 0.50 s 

As 3s'^3p^3d^°4s^4p'^4d° 0.50 0.50 0.45 p 



enables one to distinguish two contributions to the self- 
energy: a screened exchange part {'Esx) arising from the 
poles of Go, and a dynamical Coulomb interaction be- 
tween an electron and the hole-like charge distribution 
around it (Sch) stemming from the poles of Wq. The 
former tends to increase the quasiparticle energy, after 
the bare exchange is excluded. The latter tends to de- 
crease the quasiparticle energy. In particular, the dynam- 
ical contribution T,ch is highly sensitive to the number of 
bands n included in the calculation of the Green's func- 
tion in Eq. (O . The final quasiparticle energy of a state 
fnk is given by {Vxc is the LDA exchange-correlation po- 
tential) 

Kl = + (^nkl S(r, r'; E^) - F.c(r)|^„k)- (3) 



III. RESULTS 

A semicore, non-relativistic pseudopotential was gen- 
erated for Si, using the TrouUier-Martins schemeii^ For 
Ga, Ge and As, we constructed semi-relativistic pseu- 
dopotentials using the Kerker schemeiiSi These choices 
resulted in stable, transferable pseudopotentials in the 
Kleinman-Bylander form, without ghost statesii^ A sum- 
mary of atomic parameters is presented in Table A 
good expansion of electronic wavefunctions in a plane- 
wave basis was obtained using a cutoff energy of 700 Ry 
(600 Ry for Si), and the first Brillouin zone was sampled 
using a 4x4x4 Monkhorst-Pack grid.^^ These numerical 
parameters ensure convergence in LDA energy eigenval- 
ues to 0.01 eV or better. The Ceperley- Alder exchange- 
correlation potential is used.-li' For the lattice parameter, 
we used the experimental values: 5.43 A, 5.65 A, and 5.66 
A for Si, Ge, and GaAs, respectively. The polarizability 
was expanded in a plane-wave basis with an energy cut- 
off of 45 Ry (50 Ry for Si) and numerically inverted for 
the calculation of the screened Coulomb interaction Wq. 
Numerical precision in the calculation of the self-energy 
is 0.05 eV or better. 

Table shows some of the energy gaps obtained 
in the present approach for Si, compared with previ- 



3 



TABLE II: Bands Gaps of Si. All quantities in eV. 
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100 



150 



200 



Egap 


I — 1 


1 — A 


LDA present work 






0.46 


2.52 


0.60 


GW valence pseudopotential + CPP 






Shirley et al. 1.13 


3.28 


1.31 


GW present work 






1.04 


3.24 


1.18 


GW all-electron 






Hamada et al} 1.01 


3.30 


1.14 


Kotani and van Schilfgaarde'^ 0.89 


3.12 




Ku and Eguiluz'* 0.85 


3.12 




Experiment*^ 1.17 


3.35 


1.3 



" Reference 9. 
Reference 5. 
Reference 6. 

non-self-consistent results from Reference 7. 
Reference 23. 



ous valence-only pseudopotential and all-electron calcu- 
lations. Overall agreement between the present results 
(which explicitly include the effect of the semicore states) 
and experimental measurements is at the level of 0.1 
eV, and discrepancies between our results and previous 
pseudopotential-based calculations of Ref. 9 are equally 
small. Recent all-electron calculations carried out at the 
same level of the GW approximation, however, system- 
atically underestimate the minimum gap and the direct 
r — r gaps&i We find that the convergence of the self- 
energy with respect to the number of unoccupied bands 
included in Gq in Eq. ^ is an important factor. In 
Figure n we show the behavior of the calculated energy 
gap as function of the number of unoccupied bands, ric, 
included in Go- Convergence is typically very slow, and 
well-converged results require Uc > 120. Other energy 
transitions (T — X and F — F) show similar behavior 
and also approach the converged value from below. In 
contrast, the results of Ref. 7 were obtained with only 
ric = 24 and are closer to our results at approximately the 
same value of ric than to the converged results, as shown 
in Figure n This fact points to a lack of numerical con- 
vergence in the evaluation of the self-energy corrections 
in Ref. 7. 

The energy gaps obtained for Ge are presented in Ta- 
ble Spin-orbit interactions are included as first-order 
perturbations.^® A common feature of LDA-based calcu- 
lations is the overlapping of the valence and conduction 
bands at the F point. The inclusion of 3d electrons moves 
the conduction bands further down resulting in a sizable 
negative direct gapiS*2*24iiiS This feature is verified in our 
LDA calculation. As shown in Table ITTTl self-energy cor- 
rections are responsible for an opening of the gap and the 
correct positioning of the minimum, indirect gap between 
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FIG. 2: Similar to Fig. 0for Ge. Numerical values for the 
gaps r — r (upper panel) , T — L (middle panel) , and T — X 
(lower panel) are shown in open triangles. The full line is a 
guide to the eye. Results obtained by Ku and Eguilua- are 
shown in black diamonds. 



TABLE III: Bands Gaps of Ge. All quantities in eV. 



F" - 




F" -X" 


LDA present work 

-0.04 


-0.26 


0.56 


GW valence pseudopotential -I- CPP 
Shirley et al. " 0.73 


0.85 


1.09 


GW present work 

0.65 


0.85 


0.98 


GW all-electron 

Kotani and van Schilfgaarde'' 0.47 
Ku and Eguiluz'^ 0.51 


0.79 
1.11 


0.49 


Experiment'* 0.74 


0.90 


1.3 



" Reference 9. 
Reference 6 after inclusion of spin-orbit effects, 
non-self-consistent results from Reference 7. 
Reference 23. 
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TABLE IV: Bands Gaps of GaAs. All quantities in eV. 
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1 — 1 
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LDA present work 






0.13 


0.70 
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GW valence pseudopotential + GPP 






bnirley et al. iAZ 


1.75 


1.95 


GW present work 






1.38 


1.65 


1.83 


GW all-electron 






Kotani and van Schilfgaarde'' 1.20 


1.40 


1.46 


Experiment'^ 1.52 


1.815 


1.98 



" Reference 9. 

^ Reference 6 after inclusion of spin-orbit effects. 
" References 23,24. 



points r and L in the Brillouin zone. Our GW results 
compare well with experiment, although the T — X gap 
still shows a large discrepancy. On the other hand, there 
are significant differences between the pseudopotential 
and all-electron gaps at this level of the GW approxima- 
tion, in particular regarding the T — X gap. As found in 
previous studies, the self-energy corrections are needed 
to give the correct band topology. 

Figure 121 shows direct and indirect energy gaps in Ge 
as function of the number of unoccupied bands included 
in the calculation of Gq. Whereas the indirect gaps 
T — L and T — X approach the converged value from 
below, the direct gap F — F approaches it from above. 
This particular convergence behavior arises from the fact 
that we are plotting differences of quasi-particle energies: 
taken individually, all quasi-particle energies E"^^ con- 
verge monotonically from above, reflecting the attractive 
nature of the Coulomb-hole Y^ch term.^-^ Additionaly, we 
note that in Ref. 7 the direct gap at F is strongly over- 
estimated and the indirect gap F — X is underestimated 
by « 0.8 eV. As expected, the same pattern of overesti- 
mation/underestimation is evident in Figure |21 when we 
reduce the number of unoccupied bands from 170 to 24, 
which was the value used in Ref. 7. 

GaAs shows behavior similar to Ge. Table Hvl summa- 
rizes our results. Agreement with experimental data is 
now within 0.15 eV, and all-electron calculations again 
underestimate the energy gaps. Regarding the conver- 
gence with ric, Figure |3| shows that the indirect gaps 
converges slowly from below. On the other hand, the 
direct F — F gap converges more quickly. 

In all systems studied, we applied the GPP model^ in 
two ways: (1) using only the valence charge density in 
the f-sum rule, and (2) using a total charge density from 
both the valence and semicore bands. Method (1) is phys- 
ically more realistic, and predicts an energy dependence 
of the inverse dielectric function that is consistent with 
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FIG. 3: Similar to Fig. I^jfor GaAs. The thick horizontal line 
represents results obtained by Kotani and van Schilfgaarde- 
(the number of bands used was not reported in this reference). 



the RPA. Method (2) implies that semicore electrons are 
able to screen electric fields as efficiently as valence elec- 
trons, which is not physical. Nevertheless, we find agree- 
ment between the two schemes to better than 0.1 eV in 
the converged energy gap for all systems studied. This 
is a consequence of cancellation of errors: in method (2), 
the plasma frequency is overestimated, and therefore the 
contributions T,sx and Sch are overestimated in absolute 
value. Since they have opposite sign, the final quasipar- 
ticle energies are weakly affected. We also observe that 
method (2) shows slower convergence of energy gaps with 
respect to ric, which is expected since the Sch term is en- 
hanced. All results presented in this article were obtained 
using the physically more appropriate method (1). 



IV. DISCUSSION 

As shown in Figures fll3l energy gaps in the GW ap- 
proximation have a significant dependence on the num- 
ber of unoccupied bands that is slow to converge in many 
cases. This fact was already reported in calculations us- 
ing standard pseudopotential techniques. ^'^^ We can ana- 
lyze the physics of this convergence by examining the GW 
approximation in the static limit, the so-called COHSEX 
approximationiSii^ Bands gaps can be calculated more 
easily within the COHSEX approximation and show con- 
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vergence behavior similar to the full dynamic calcula- 
tion. Under this approximation, the Coulomb-liole term 

I 



T.coH{vy) = -&{v~v')Wpoi{v,v') 



where the second equality follows from completeness of 
the basis of eigenvectors. In actual calculations, this sum 
over bands n is always truncated, and the equality is 

I 



where Wgg' (q) are the coefficients in the plane- wave ex- 
pansion of the polarization potential^*^ and we define 

X^'"(k, q) = (nk - qle-^^'i+'^^-nrnk) . (6) 

Physically, the summation in Eq. describes vir- 
tual transitions produced when the quasiparticle induces 
a charge fluctuation around itself. Y,coh is the energy as- 
sociated to the interaction between the quasiparticle and 
the induced charge fluctuation. The matrix elements, Eq. 
©, decay slowly as the energy difference between bands 
m and n increases. Monitoring the convergence of these 
matrix elements provides a good estimate of the relative 
error in the T,coh, but gauging the absolute convergence 
of the self-energy requires knowing Wpoi, which depends 
on the physical system. 

Although we have not directly investigated the impor- 
tance of imposing self-consistency in the calculation of 
the self-energy, this is an unsetled issue and we address 
it briefly below. Von Barth and Holm have investigated 
the effect of self-consistency in the electron gas;^ and 
concluded that restricted self-consistency has small but 
significant effect on the full bandwidth and in the satel- 
lite structure of the electron gas. On the other hand, 
full self-consistency gives a poor description of the satel- 
lite structure and the bandwidth is drastically increased. 
Inclusion of vertex corrections are expected to recover 
the good, non-self-consistent results, but calculation of 
vertex corrections is not a simple task even for the elec- 
tron gas system. ^•^-'^ Self-consistency has been recently 
applied to real materials (iSS and the valence bandwidth 
is also shown to increase when self-consistency is im- 
posed. It appears that one must therefore include self- 
consistency and vertex corrections together in order to 
obtain a meaningful picture. Inclusion of vertex correc- 



has a simple form in terms of the polarization potential, 
Wpoi[r, r') = W[r, r', u^Q)-V{t- r'), 



- 2E^"k(rXk(r')M^po;(r,r') , (4) 

nk 

I 

violated. The Coulomb-hole energy evaluated at a given 
electronic state is calculated according to the expression 



(5) 



tions and self-consistency is tangential for the purpose of 
our work, which is to compare pseudopotential-based and 
all-electron GW calculations and understand the role of 
semicore electrons. 



V. CONCLUSION 

We conclude that for the systems considered, the 
valence-only pseudopotential method does not suffer 
from large errors from the neglect of core states, as 
claimed in some all-electron calculations.^'^ While semi- 
core effects are negligible in Si, they are important if one 
aims at good quantitative agreement with experiment 
in Ge and GaAsiSi^ii^ii^ However, discrepancies between 
pseudopotential and all-electron based GW calculations 
reported in recent worka^ may be explained by a lack 
of numerical convergence in the latter. Specifically, the 
self-energy calculated within the GW method has slow 
convergence with respect to the number of energy bands 
included in the calculation of the Green's function, as is 
demonstrated in this work and has been pointed out in 
the pastiSiSS This convergence behavior is present in the 
static limit to GW, the COHSEX approximation, and 
can be analyzed by comparing the COHSEX Coulomb- 
hole energy, Eq. obtained with and without explicit 
summation over energy bands. 



Acknowledgments 

We thank Andrew Canning and David Clatterbuck 
for fruitful discussions regarding the FLAPW method. 
This work was supported by National Science Founda- 
tion Grant No. DMROO-87088 and by the Director, Office 



(mk|Scoff|mk) = - T4-GG'(q)5][Mr(k,q)]*>('c^'7(k,q) 

GG'q n 

I 



6 



of Science, OfRce of Basic Energy Sciences, Division of 
Materials Sciences and Engineering, U.S. Department of 
Energy under Contract No. DE-AC03-76SF00098. Com- 
putational resources have been provided by NSF at the 



National Partnership for Advanced Computational In- 
frastructure, and DOE at the National Energy Research 
Scientific Computing Center. 



1 M.S. Hybertsen and S.G. Louie, Phys. Rev. B 34, 5390 

(1986). 

^ W.G. Aulbur, L. Jonsson, and J.W. Wilkins, Solid State 
Physics, eds. F. Scitz, D. Turnbull, and H. Ehronrcich, 
Academic, New York, vol. 54, 1 (2000). 

^ F. Aryasetiawan and O. Gunnarsson, Rep. Prog. Phys. 61, 
237 (1998). 

" P.E. Blochl, Phys. Rev B 50, 17953 (1994). 

® N. Hamada, M. Hwang, and A.J. Freeman, Phys. Rev. B 

41, 3620 (1990). 
® T. Kotani and M. van Schilfgaarde, Sol. State. Comm. 121, 

461 (2002). 

^ W. Ku and A.G. Eguiluz, Phys. Rev. Lett. 89, 126401 
(2002). 

® S. Lcboguc, B. Arnaud, M. Alouani, and P.E. Blochl, Phys. 

Rev B 67, 155208 (2003). 
^ E.L. Shirley, X. Zhu, and S.G. Louie, Phys. Rev. B 56, 

6648 (1997). 

R.M. Dreizler and E.K.U Gross, Density functional theory: 
an approach to the quantum many-body problem (Springer- 
Verlag, Berlin, 1990). 

J. P. Pordcw and A. Zungor, Phys. Rev. B 23, 5048 (1981). 
N. TrouUier and J.L. Martins, Phys. Rev. B 43, 1993 
(1990). 

13 G.P. Kerker, J. Phys. C 13, L189 (1980). 

" M. Rohlfing, P. Kriiger, and J. Pollmann, Phys. Rev. B 



57, 6485 (1998). 
1^ L. Hedin and S. Lundqvist, Solid State Physics, eds. F. 
Scitz, D. Turnbull, and H. Ehrenreich, vol. 23, 1, Aca- 
demic, New York (1969). 

L. Kloinman and D. M. Bylander, Phys. Rev. Lett. 48, 
1425 (1982). 

H.J. Monkhorst and J.D. Pack, Phys. Rev B 13, 5188 
(1976). 

1* M.S. Hybertsen and S.G. Louie, Phys. Rev. B 34, 2920 
(1986). 

1® F. Aryasetiawan and O. Gunnarsson, Phys. Rev. B 54, 
17564 (1996). 

^° A. Fleszar and W. Hanke, Phys. Rev. 56, 10228 (1997). 
^1 U. von Barth and B. Holm, Phys. Rev. B 54, 8411 (1996); 

Phys. Rev. B 55, 10120 (1997); B. Holm and U. von Barth, 

Phys. Rev. B 57, 2108 (1998). 

W-D. Schone and A.G. Eguiluz, Phys. Rev. Lett. 81, 1662 
(1998). 

^3 Landolt-Bornstein, Numerical Data and Functional Rela- 
tionships in Science and Technology, eds. K.-H. Hellwege, 
O. Madelung, M. Schulz, and H. Weiss, New Series, vol 
III, Pt. 22a, Springer- Vorlag, Now York (1987). 
P. Lautenschlager, M. Garriga, S. Logothetidis, and M. 
Cardona, Phys. Rev. B 35, 9174 (1987). 



